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EHRHART POLYNOMIAL FOR LATTICE SQUARES, CUBES AND 

HYPERCUBES 

EUGEN J. IONASCU 


Abstract. In this paper we are constructing integer lattice squares, cubes or hypercubes in R d 
with d £ {2,3,4}. We find a complete description of their Ehrhart polynomial. We characterize 
all the integer squares in R 4 , in terms of two Pythagorean quadruple representations of the form 
a 2 + b 2 + c 2 = d 2 , and then prove a parametrization in terms of two quaternions of all such squares. 
We introduce the sequence of almost perfect squares in dimension n. In dimension two, this is very 
close to the sequence A194154 (in OEIS). 


1. INTRODUCTION 

Eugene Ehrhart ([8], [9]) proved that given a d-dimensional compact simplicial complex in M n 
(1 < d < n), denoted here generically by V, whose vertices are in the lattice Z n , there exists a 
polynomial L(V,t) E Q[t] of degree d, associated with V, satisfying 

L(V,t) = the cardinality of {tV} n Z n , t € N. 

It is very interesting that one can say more about three of the coefficients of L(V, t): 

(1) L(V,t) = Vol(V)t d + ^Vol(dV)t d ~ 1 + ... + X (V), 

where Vol(V) is the usual volume of V normalized with respect to the sublattice containing V, 
Vol(dV) is the surface area of V normalized with respect to the sublattice on each face of V and 
x(P) is the Euler characteristic of V (in the sense of polytopal complexes, so for a convex polytope 
it is equal to one). In general, the other coefficients in dTJ) may have complicated expressions in 
terms of the vertices of V. 

It is also known that the number of points in the interior of tV is given by (—1 ) d L(V, — t ). For 
a polytope that is a cross product of two polytopes, the Ehrhart polynomial is the product of the 
corresponding smaller degree Ehrhart polynomials. In this article, we are studying this polynomial 
for lattice squares, cubes, and hypercubes in M n . Such objects have been constructed and counted 
in several works (see m, m, m, oso- 
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Figure 1. Square with a=5 and b=2 

2. Squares 

2.1. Squares in two dimensions. How general can a lattice square in two dimensions look like? 
We may assume that one of the vertices is at the origin. Then the other vertices, in counterclockwise 
order, are given by A(a, b ), B(a — b, a + b) and C(— 6, a) where a and b are integers, not both zero, 
and such that gcd(a,6) = 1. Using the facts mentioned in the Introduction, the polynomial in (|TJ) 
is simply E\j(t) = (a 2 + b 2 )t 2 + 2t + 1 and of course, EL (£) = ( a 2 + b 2 )t 2 — 2£ + 1, £ € N. 

In Figurc[I] we see the 28 points inside of the square OABC , with A{ 5, 2), B( 3,8) and C(— 2, 5). 
One interesting problem here is to determine the sequence of possible lattice points in the interior 
of such a square. The first hundred terms of this sequence, listed in increasing order, are included 
in the table below: 
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220 

221 

225 
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232 

240 

241 

249 

256 

257 

264 

265 

268 
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276 
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288 

289 

292 

293 

297 

301 

304 

305 

312 

313 

316 

324 

325 

328 

336 

337 

345 

348 


Let us call this sequence the almost perfect squares sequence. In Figure [2] we see some of the 
squares that show that the above numbers are indeed in the sequence. As a curiosity, 2015 is 
not in this sequence but 2016 is. This sequence is very close to A194154 in the OEIS (On-Line 
Encyclopedia of Integer Sequences), but 20 is the first term that is not in our sequence. 
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Figure 2. Various squares 




Let us give a reformulation of the exact form of the Ehrhart polynomial above in conjunction 
with the reciprocity property (about interior points), which can by shown independently using 
Pick’s Theorem. 


Proposition 2.1. Given a, b with gcd(a, b) = 1 and t, € N, then 

j]{(x, y) € Z 2 : ax + by , ay — bx € [1, t(a 2 + b 2 ) — 1]} = (a 2 + b 2 )t 2 — 2t + 1. 


Proof. Counting the points inside the square OABC , with 0(0,0), A(a,b), B(a — b,a + b ) 
and C (— b, a), is equivalent to count the solutions (x, y) of the system (x, y) = a(ta, tb ) + /3(—tb, ta) 
with a, /3 € (0, 1 ). Since a = (x, y) ■ (a, b)/t(a 2 + b 2 ) = and j3 = (x, y ) • (-6, a)/t(a 2 + b 2 ) = 

t(a^+y 2 ) we see cons t ra i n t s on a and j3 is equivalent to ax + by, ay — bx € [1, t(a 2 + b 2 ) — 1]. 

Hence, the result follows from the Ehrhart polynomial expression. ■ 

So, every term of the sequence of almost perfect squares is the answer to a counting as in 
Proposition 12.11 which in particular gives the justification of one of our earlier claims: 


U{(x, y) € Z 2 : 44x + 9 y, 44 y - 9x £ [1, 2016]} = 2016. 

2.2. Squares in three dimensions. In M 3 , we can obtain a lattice square by taking two orthog¬ 
onal vectors with integer coordinates of the same length: u = (a, b, c), v = ( a',b',c ') such that 
a 2 + b 2 + c 2 = a' 2 + b' 2 + c' 2 = l and aa' + bb' + cc' = 0. In [12], two such vectors are call twin 
vectors. The number of such twin vectors having a given length is calculated in [12] • In [T9] , this 
concept is generalized to m dimensions and called m-icube (a set of m vectors in Z n , mutually 
orthogonal and of the same length). Since we are interested in generalizing these ideas to any 
dimension, we are going to use this terminology also. 
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Although it does not make much of a difference it is natural to assume that the square is 
irreducible , i.e., gcd (a,b, c, a',b', d) = 1. If we set d = gcd (a,b,c), d' = gcd(a' ,b' , d) and D = 
gcd (be' — b'c , ad — a'c , ab' — b'a) , we have the following more general formula, than in two dimensions, 
for the Eu{t). 

Theorem 2.2. The Ehrhart polynomial of a lattice square embedded into R 3 , described above 
and with the notation introduced is given by 

(2) E\j(t) = Dt 2 + (d + d')t + 1. 

Proof. The fundamental domain of the lattice containing the triangle has a volume equal to 
\/n\ + 7 I 2 + n 2 where n = [n \, 712 , ^ 13 ] is a vector normal to the plane containing the triangle and 
such that gcd(ni, n 2 , 77 , 3 ) = 1 (see [lj and [18]). Such a vector is given clearly by the cross-product of 
u and v: n = u x v = j){bc' — b'c , — (ad — a'c),ab' — b'a). Let us observe that because of Lagrange’s 
Identity £ 2 = (a 2 + b 2 + c 2 )(a' 2 + b' 2 + c' 2 ) = (aa'+ bb'+ cc') 2 + (bd— b'c) 2 + (ad— a'c) 2 + (ab' — b'a) 2 , 
we conclude that |n| = jj. Hence, we see that the first coefficient of the Ehrhart polynomial is 
Vf/(i) = D. Similar argument goes for the second coefficient. ■ 

Standard Examples: Let us analyze some examples which are not xy, xz or yz-plane examples. 
First, let us take u = (3, —3,0) and v = (1,1,4); observe that |ii | 2 = \v\ 2 = 18 and u ■ v = 0. In this 
case n = 6 (— 2, —2, 1), and the equation of the plane is 2x + 2y — z = 0 and £h(f) = 6 t 2 + 4t + 1. 
We know that En(—1) = 3 represents the number of lattice points in the interior of the square. 
We notice that 3 (and Eu{— 2) = 17) was not an almost perfect number. Hence we may want to 
generalize that sequence to almost perfect squares in dimension n. 

What is interesting is that in the plane 2x + 2y — z = 0 we have another square: u = (2, —1, 2) 
and v = (—1,2,2). Also, the square above can be written in terms of this one basically like in 
two dimensions: u = u — v and v = u + v. For this square, n = 3(—2, —2,1) and so the Ehrhart 
polynomial is E{u V y(t) = 3 1 2 + 2t + 1. We see that Es^^{— 1) = 2 which is yet new for the almost 
perfect square sequence compared with dimension two. 

If we take the twin vectors u = (6, —2,3) and v = (—2,3, 6) we get n = 7(—3, —6, 2) and so the 
Ehrhart polynomial is E(t) = 7f 2 + 2t + 1. There seems to be very few numbers that are not perfect 
squares in 3D: 7, 14, 23,... 

For an example of an irreducible square, for which d and d' in Theorem 12.21 are both greater than 
one, we refer to u = 5(8,12, 9) and v = 17(0, —3,4). Its polynomial is E(t) = (5t + 1)(17 1 + 1). 

How does one construct such squares in M 3 ? We observed in the proof of Theorem 12.21 that every 
such square is contained in a plane whose normal n = {n\,n 2 ,n^) satisfies 

(3) t 2 = n\ + + n\. 
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A solution of this type of Diphantine equation is usually referred in the literature as a Pythagorean 
quadruple. For the number of primitive solutions of 0 in terms of £, we refer the reader to a recent 
paper of Werner Hiirlimann ([H]) but also [3]. There are at least two other terms used for these 
quadruples of integers: cuboids (M) and Lorenz quadruples (see m)- We include in the next 
tables all primitive solutions for odd values £ € {1,3, ..19}: 


£ 

a,b,c , gcd (a,b,c) = 1, c even 

£ 

a,b,c , gcd(a, b,c) = 1 , c even 

1 

[1,0,0] 

11 

[9, 2, 6], [7, 6, 6] 

3 

[1,2,2] 

13 

[5, 12, 0], [3, 4, 12] 

5 

[3,4,0] 

15 

[5, 14, 2], [11, 2, 10] 

7 

[3,6,2] 

17 

[15, 0, 8], [1, 12, 12], [9, 12, 8] 

9 

[1, 4, 8], [7, 4, 4] 

19 

[1, 6, 18], [17, 6, 6], [15, 6, 10] 


So, having a Pythagorean quadruple, how can we construct a related twin set of vectors? We 
include a partial answer to this question at this point, but we will include the complete solution 
later in this paper. 

Theorem 2.3. Suppose thatu' = (o, b, c) satisfies n\a + n?b + n^c = 0 , where n\ + n 2 + n 2 = t 2 
with all variables involved being integers. Then there exist v = ( a',b',c') such that tv! and v define 
a lattice square in the plane of normal n = (m, 712 , 71 , 3 ). 

Proof. We define v to be the cross-product of u' and n. Clearly, v is perpendicular to n and 
so it is in the right plane. Since v! and n are perpendicular, |u| = |t/'11| = £\u'\. ■ 

It is natural to look for the “smallest” square in the plane nia + 77 - 2 ^ + R 3 c = 0. This leads 
us to the shortest vector problem (SVP), i.e., finding a non-zero vector in a lattice of minimum 
norm. An interesting problem, at this point, is to characterize all the values l so that Vt appears 
as side-lengths for an lattice embedded square in R m , m > 2. If we denote this set by £, from 
the two dimensions construction, we see that C is invariant under multiplication with numbers 
which are sums of two squares. For m = 3, it is also clear all the numbers of the form 4 fc (8s — 1) 
are not in C, since these are not representable as sums of three squares (Legendre’s three-square 
theorem). Also, 3, 11 or 19 (see the sequence A223732 in OEIS) are not in C, since they have only 
one representation as sum of three squares and these representations contain only odd numbers. 
Hence the sum aa' + bb' + cd is also odd and so, it cannot be zero. We will describe the set L at 
the end of the subsection. We will see that in M 4 and beyond, C = N. 

We have the following parametrization for those squares whose side-lengths are natural numbers 
of the form x 2 + y 2 + z 2 + t 2 (see [T7] and [23]): 


(4) u := (2 ty + 2 zx, 2 tz — 2 yx, t 2 — z 2 — y 2 + x 2 ), v : = (2 zy — 2 tx, z 2 — t 2 + x 2 — y 2 , 2tz + 2 yx), 
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having normal vector n = (— x 2 + t 2 — y 2 + z 2 , —2(tx + zy), 2{ty — zx)), |n| = x 2 + y 2 + z 2 + t 2 . 
Next, we show that, as in m and [241 . the following similar parametrization for the solutions of 
@ takes place. The proof is essentially the same as in [24] but we include it for the convenience of 
the reader. 

Theorem 2.4. f[24j ) Every primitive solution of after a permutation of variables and 
change of signs, is given by n\ = 2 {zy — tx), n 2 = 2 (tz + yx) and n 3 = z 2 — t 2 + x 2 — y 2 , and 
l = x 2 + y 2 + z 2 + t 2 for some integers x, y, z and t. 

Proof. As usual, a primitive solution of © is one for which gcd(ni, n 2 , n 3 ) = 1. We must 
have i odd since otherwise gcd(ni, 71 - 2 , n 3 ) > 2. In this case, one of the n* must be odd and the 
other two even. Without loss of generality, let us assume that n\ and ri 2 are even. Then, we have 

(£ - n 3 )(i + n 3 ) = n\ + n 2 . 

We know that i + n 3 = 2 a and l — n 3 = 2(3 are both even, and then the above equality becomes 

a/3 = (m/2) 2 + (n 2 /2) 2 = A 2 + B 2 = (A + Bi)(A - Bi ), 
where A = n\/2 and B = n 2 /2. A Gaussian prime of the form p = 4k + 3 which divides a/3, must 
divide both A and B. So, it cannot divide both a and /3 because it then divides i = (a + (3)/2 
and n 3 = (a — /3)/2, which implies that gcd(ni, ?r 2 , n 3 ) > p. Hence, a and f3 are both sums of two 
squares. If we let gcd(A+ /.£>, a) = t + iy and gcd(A + iH,/3) = x + iz. Let us make the observation 
that these equalities are defined up to a unit, i.e., ±1 or ±i. 

Since a is real, t — iy divides a and so a = (t 2 + y 2 )a'. But every prime factor of a, appears as 
a factor of either A Bi or A — Bi. Taking into account the multiplicities we see that a divides 
t 2 + y 2 and so a' = 1. Similarly, we have j3 = x 2 + z 2 . Also, t + iy divides A + iB and also x + iz 
divides A + iB. Hence A + iB = (t + iy)(x + iz)k and from here A 2 + B 2 = (t 2 + y 2 )(x 2 + z 2 )\k\ 2 
which together with what we have shown earlier forces |£;| = 1. So, by changing x and y we can 
assume that k = 1. Hence we have A = tx — zy and B = xy + tz. Then l = a + /3 = x 2 + y 2 + z 2 +t 2 
and n 3 = a — /3 = t 2 + y 2 — x 2 — z 2 . ■ 

Theorem 12.41 allows us to characterize C. The result is not new, as we found recently, it is 
contained in [12j. We include a proof of it, based on our development of the ideas. 

Theorem 2.5. The set of all l so that VI is the side-lengths for an embedded square in Z 3 is 
the set of positive integers which are sums of two squares. 

Proof. In one direction, i.e., C contains the set of all positive integers which are sums of two 
squares, we can observe that the lengths for squares in two dimensions are of the form \/a 2 + b 2 . 
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For the other direction, if we start with an arbitrary square, 5, in Z 3 , let us suppose the side is 
\ft. As in the proof of Theorem 12.31 this square is contained in a plane of normal n = (n-i, U 2 > R 3 )> 
with n\ + n 2 + n 2 = l 2 - By Theorem 12.41 we can find x, y, z and t, such that n\ = 2 zy — 2 tx, 
ri 2 = 2 tz + 2 yx and n% = z 2 — t 2 + x 2 — y 2 . Then the parametrization ((4]) gives a square in the same 
plane as S, whose sides are £. We can then write the vector of this square in the basis given by two 
of the vectors in S: w = au + /3v for some a and (3 € Q. Then, taking norms we get f' 2 = a 2 £ + /3 2 £ 
so t = a 2 + /3 2 . This shows that £ = A 2 + B 2 for some A, B £ Z. ■ 

Theorem 12.31 and Theorem 12.51 imply the next consequence which one can also prove elementary. 

Corollary 2.6. Given a point P = ( x,y,z ) € Z 3 in the plane of equation n\x + n^y + n^z = 0 
with n\ + n 2 + n 2 = l 2 for some < £ N, then the number x 2 + y 2 + z 2 is actually a sum of two 
squares. 

2.3. Squares in R 4 . Perhaps, one of the simplest ways to construct squares in Z 4 is to take two 
squares in two dimensions and “add” them together, each on its own dimensions. In other words, 
the two vectors u and v that define the square, as we have seen before, are of the form u = (a, b, c, d) 
and v = (—6, a, —d, c ). We have clearly u ■ v = 0 and ||u|| 2 = ||u|| 2 = a 2 + b 2 + c 2 + d 2 . This simple 
example allows us to answer the question about the possible side-lengths of such a squares. By 
Lagrange’s Four Square Theorem, we see that every natural number is a possible side-length of a 
square. Are these squares the most general situation that one can expect? We observe that u and 
v are essentially the first two rows of the following pseudo-orthogonal matrix (every two rows are 
mutually orthogonal and have the same norm): 

abed 
n _ —b a d —c 

Ua ’ b ’ C ’ d - - c -d a b ' 

—d c —b a 

It is not difficult to see that the product of matrices like these ( O t O = tl ) is of the same form. 
This allows one to define a certain multiplication on vectors in R 4 which is exactly the quaternion 
multiplication that we will be using latter. 

Example 1. Let us compute the Ehrhart polynomial for a particular case like this that which 
gives a new value for its sides. We have u = (2,1,1,1) and v = (—1,2, —1,1). Using the same idea 
as in the three dimensions we are computing two normal vectors that define the orthogonal space 
of the plane generated by u and v. For a generic point P = ( x,y,z,t ) € R 4 , the two equations, 
computed with the help of the cross-product in the three dimensions, are 


3x — y — 5z = 0 and 2y + z — 3f = 0. 
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We observe that both u and v satisfy these equations. Since these two equations can be solved for 
y and t, we can find two integer vectors which generate the lattice Z 4 intersected with the plane 
defined by u and v: y = 3.x — 5 z and t = 2x — 3 z which gives the two vectors a = (1, 3,0, 2) and 
P = (0, —5,1, —3). We know then that these two vectors form a fundamental domain, so we need 
to compute the area of the parallelogram generated by these two vectors: A = |ck| |/9| sin('y) where 
cosy = a ■ /5/|ck||/3| = —3/t/IO. This gives A = 7 and so, the Ehrhart polynomial is 

= t 2 + 2t + 1 = (t + l) 2 . 

For a set of vectors S in M 4 we denote as usual by 5^ the set of all vectors x € M 4 perpendicular 
to every vector v in S, i.e., 

S 1 - = {(xi,X 2 ,x 3 ,x 4 ) G M 4 |xiWi +X 2 V 2 + X 3 V 3 + X 4 V 4 = 0 for all v = (vi, v 2 , v 3 , V4), v G 5}. 

Theorem 2.7. ( i ) Given integer vectors u = (111,112,113,114) and v = (ui, V2, V3, V4) such that 

0 < t = u\ + ul + u\ + u\ = v\ + v\ + v\ + v\, UlVl + U2V2 + U3V3 + U4V4 = 0, 

then there exist an odd k G N dividing i and two vectors w\ and W2 with integer coordinates such 
that 

w\ = (0, «i - fti, a 2 - @2, «3 - P3) and w 2 = (a 3 — /3 3 , -a 2 - @2, «i + Pi,0) with 

(5) k 2 = af + 0.2 + a 3 = Pf + 02 + P3, 

ati and Pi of the same parity and u,v G {wi,W2} ± - One can permute the coordinates of u and v 
and/or change their signs in order to have w\ and W2 linearly independent. 

(ii) If gcd(ai, Pi, a2, P 2 , P3) = 1 then the volume of the fundamental domain of the minimal 

lattice containing u and v is equal to k. 

Proof. Using Lagrange’s Identity, we get 

44 4 

t 2 = E = E UiVi ) 2 + (■ Uiv i - u 3 v if- 

i=l i= 1 i= 1 l<«<i<4 

Since uiv 1 + U2V2 + U3V3 + U4V4 = 0, this implies that 

( 6 ) £ 2 = ^2 (uiVj - Ujv/ 2 . 

l<i<j<4 

If we denote by A tJ = (—1 (uiVj — Ujvf), it is not difficult to check that 
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( 7 ) A12A34 — A13A24 + A14A23 — 0 . 

This helps us change dU) into 


( 8 ) 


ft = (A12 ± A34) 2 + (A13 =F A24) 2 + (A14 ± A23) 2 . 


One can check that 


ui( 0 ) + U 2 A 34 + U 3 A 24 + W 4 A 23 = 0 
/„N ( Ul( 0 ) + U 2 A 34 + U 3 A 24 + U 4 A 23 = 0 

U 1 A 23 + U 2 A 13 + U 3 A 12 + 114 ( 0 ) = 0 
^iA 23 + v 2 A V3 + U 3 A 12 + u 4 (0) = 0. 

Let us observe that if 2 divides i, then all A ij are divisible by 2 , because all Ui (vi) are all even 
or all odd. In order to prove the claim, we simplify ([S]) by the greatest power of 2 possible and 
take a\ = (Ai 2 + A 34 )//, a 2 = (—A13 + A 2 4 )/fa 0:3 = (Ai 4 + A 2 3)/i', fa = (Ai 2 — A 34 )/£', 
fa = (-A13 - A24 )/fa and fa = (A u - A 23 )/£', with £' so that gcd(«i, fa, a 2 , fa, a 3 , fa) = 1 . 

It is easy to check that if «3 7 ^ fa then the two vectors w\ and w 2 are linearly independent and 
so the orthogonal space is two-dimensional. By permuting coordinates and/or changing 

their signs, we can insure that 03 — fa = 2 A 2 3 /£' / 0 , since by (| 6 |) and the assumption that £ > 0 , 
we see that not all of the A l? are equal to zero. 

(ii) Let us denote by C the minimal lattice which contains the given square {u, u} and by V its 
volume. As we observed in part (i) this lattice is the same as Z 4 n if a3 / fa- Let us 

observe that the two vectors U := (0, A 12 , — A 13 , A 14 ) and V := (A 14 , — A 24 , A 34 , 0) are in C. This 
shows that 2 U/£' = (0, ot\ + fa,a 2 + fa, 013 + fa) and 2V/Z' = ( 0:3 + fa, —a 2 + fa, ot\ — fa, 0) are in 
C also. The area of the parallelogram determined by these two vectors is given by the square root 
of the Gramian determinant 

(«i + fa ) 2 + («2 + fa ) 2 + («3 + fa ) 2 («1 + fa)(fa - a 2 ) + (a 2 + fa)(ai - fa) 

(ai + fa) {fa — a 2 ) + (a 2 + fa){cti — fa) ( 0:3 + fa ) 2 + ( a 2 — fa ) 2 + (aq — fa ) 2 

_ 2 {k 2 + a.i fa + ot 2 fa + Oi 3 fa) 2 {oi\fa — faa 2 ) 

2 {a\fa — facx 2 ) 2 {k 2 + oi 3 fa — a 2 fa — oiifa) 

This means that the vectors U/£' and Vj fa still in C, form a parallelogram of an area, which is the 

square root of 


(1/ 4 ) [( k 2 + 013 fa) 2 — ( a 2 fa + a±fa) 2 — (aifa — faa 2 ) 2 ] 
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(1/ 4 )[k A + 2 k~otzfa + ot 2 ( 3 2 — ( k 2 — a 2 )(k 2 — / 3 f)] — /c 2 ((«3 + fa)/2 ) 2 . 

It follows that V divides fc |«3 + /3s|/2. In a similar way we can show that 1/ divides =t |/2 for 
i = 1,2,3. Given the assumption that gcd(ai, fa, 02 , fa, « 3 , fa) = 1, we conclude that V divides k. 

In order to conclude that V = k, let us look at the Gram determinant of the vectors w\/{013 — fa) 
and 102/(03 ~ fa), which if added to £ extends the lattice and gives a basis for it. Since its volume 
can be determined by the same method as above, a similar calculation gives V' = k/\a^ — fa\. 
This implies that V/V' is an integer, or k divides I7|ct!3 — fa\. Since we can obtain similarly that k 
divides V\a.i ± fa\, we conclude that V = k. ■ 

Example 2: As in the case of equilateral triangles (see BSI), there is a converse of Theorem l2.7l 
Let us take the first odd l for which we get two essentially different representations as in ([5]): 

ll 2 = 9 2 + 6 2 + 2 2 = 7 2 + 6 2 + 6 2 . 

We can then choose a\ = 9, /3i = 7, cx .2 = 6 , fa = 6 , 03 = 2, and fa = 6 . Then the two vectors, 
defined in Theorem l2.71 are w\ = (0,2,0, —4) and W 2 = (—4, —12,16,0). Then the space {ioi, 
is defined by the equations y — 2t = 0 and x + 3y — 4z = 0. These can be simplified to y = 2t and 
x = 4z — 6 1. Because the generic vector in {ioi, H Z 4 is u = (4 z — 6 1, 2 1, z, t ) = z( 4,0,1,0) + 
t(— 6 , 2 , 0 , 1 ) we can calculate, as before the “volume ” of the fundamental domain and obtain indeed 
11. Then we get a quadratic form QF(z, t ) = (4z — 6 1) 2 + (2t ) 2 + z 2 + 1 2 = 17 z 2 — 48 zt + 41t 2 that 
should lead to the solutions we need solving the Diophantine equation QF(z,t ) = 11 k. It turns 
out that the smallest multiple for which we have solutions is £ = 11 and then two of the vectors 
which define the square are u = (—4,8,5,4) and v = (10,2,4,1). Its Ehrhart polynomial is then 
E\j(t) = lit 2 + 2 f + 1 . We observe that we get new terms, such as 10, 94, and 266 for instance, 
for the sequence of almost perfect squares in dimension 4, compared to dimension 2. This is an 
example when all the are divisible by 11 . 

As we have seen in (fl6j). the following result gives a certain converse to Theorem 12.71 

Theorem 2 . 8 . Given k odd, and two different representations 

k 2 = a 2 + b 2 + c 2 = a' 2 + b' 2 + c' 2 , with gcd(a, b, c, a',b',c!) = 1, d > c, 

and a, a’ both odd. Then if we set A12 = , A34 = ^ 4 ^-, A13 = — A24 = ^y L , A14 = 

and A23 = the equations & and ( 0 ) are satisfied. Moreover, the two dimensional space S of 
all vectors [u, v, w, t] E Tfi , such that 


(10) 


(0 )u + A34V + A24 W + A231 = 0 
A23 U + A 13 V + A12 w + (0 )t = 0 
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contains a family of squares. 

Proof. By construction we see that © is true. Since a, a' are odd, the others must be even and 
so, all numbers defined above are integers. We have then indeed 

A 12 A 34 — A 13 A 24 + A 14 A 23 = (1/4) [a ' 2 — a 2 — (b 2 — b' 2 ) + (c /2 — c 2 )] = 0. 

One can also check that k 2 = Yli<j A 2 j. We see that the assumption A 23 > 0 insures that the 
equations (flOl) define a two dimensional space in M 4 . We can solve the equations (flOl) for t and u: 


A u v + A 24 W’ Ai 3 u + Ai 2 u> ... ^ „ 

t = - 1 -, u = - 1 -, with v,w € A. 


A 


23 


A' 


23 


If we denote a generic point P E R 4 in the plane (1101) . in terms of v and w, i.e., 


A13U + A12 W A34U + A24 w 

P(v,w) = [ --r- ,v,w, -t-J, 

-^23 ^23 

we observe that for two pairs of the parameters, (v,w) and (v',w'), we obtain A 23 = wv' — vw', 

A t , Ai 3 u' + Ai 2 w' A 13 U + A 12 W , A A, / A 
A 12 = -^-) + -X- V = Al 2 A 23 /A23, 


23 


A 


23 


and similarly A ' 13 = A 13 A 23 /A 23 , etc. 

We want show that a non-zero square the form OP(v, w)P(y',w')Q ( OQ = OP(v, w)+OP(v',w')) 
exists for some integer values of v, w, v' and w'. We need to have OP(v,w ) 2 = OP(v',w') 2 = l 
and P(v,w ) • P(v',w') = 0. 

If such a square exists, by Theorem 12.71 we may want to l of the form £ = kl' for some P € N. 

By Lagrange’s identity, calculations similar to those in the derivation of (f 6 |) , show that if 
\P(v,w)\ 2 = \P(v', w')\ 2 = kP then 


(P(», w> -P(v’, w')f = - (X; A|) (If ) :2 = k 2 (* - e (AT) 2 . 

' l iJ 

Hence, to get P(v,w ) • P(v',w') = 0, it is enough to have A ' 23 = P A 2 3 (provided that we have 
already shown that \P(v,w )\ 2 = \P{v',w ')\ 2 = kP). In other words, a square exists if there exist 
integer pairs (v,w) and ( v',w') such that 


( 11 ) 


\OP(v, w)\ 2 = | OP(v',w')\ 2 = k^^-. 

A23 


The essential object in this analysis is the quadratic form 

f A 34 u + A- 24 W \ 2 


QF(v,w) := \OP(u,v)\ = 




A' 


23 




+ v 2 + w 2 + 


( A13U + A12UI^ 2 




A' 


23 


or 
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QF(v,W ) = 


(A 34 + A 43 + A 23 )f + 2(A34A24 + Ai3Ai2)l , 'U’ + (A 24 + A 42 + ^ 23 ) w 

A 2 

^23 


We use Lagrange’s identity again, which we will write in the form 


(a 2 + /3 2 + 7 2 )(C 2 + V 2 + 9 2 ) = («C - Pv + 7#) 2 + (ay + /?() 2 + (a<9 - 7 C ) 2 + (/?# + yr?) 2 . 

Then the determinant of the form QF(v, w) (excluding its denominator) is equal to 

—A /4 = (Ajj 4 + Aj 3 + A 23 )(A 42 + A 24 + A 23 ) — (A34A24 + A13A12) 2 = 

(A12A34 — A13A24 + A 23 ) 2 + (A12A23 — A34A23) 2 + (A13A23 + A24A23) 2 => 
—A /4 = A 23 [(A23 — A14) 2 + (A12 — A34) 2 + (A13 + A24) 2 ] = A 23 (A: 2 ) 

=► A = ~( 2 kA 23 ) 2 . 

This implies that for v 3 = — (A34A24 + A13A12) and wq = A 2 3 + A 23 + A| 4 we get 

QF(v 0, wq) = k 2 ( Ag 4 + A 43 + AI3), 

and 


( 12 ) 


QF(y,w) = 


(wqv — vq w ) 2 + k 2 iu 2 A 23 
A 2 3 u;o 


So, by (HU) and (fT 2 l) we need to have solutions (v,w) and (v',w') of 


(13) {wqv — vqw ) 2 + k 2 w 2 A 23 = (wqv 1 — vqw ') 2 + k 2 w' 2 A 23 = kA 23 wo(wv' — vw'). 

In general, if the Diophantine quadratic equation x 2 + y 2 = n has a solution (x,y), then it has 
other solutions such as (y, —x ). 

Given (v, w ), let us assume that (y' , w ') gives exactly this other solution (x = wqv—vqw, y = kwA 23 ) 
when substituted in (fl3l) : 


wqv' — vo w' = kwA 23 

< 

kA 23 w' = —{wqv — vqw). 
This allows us to solve for v' and w' in order to calculate A 23 


v = 


VqW — vqWqv + k 2 A 2 3 w 


w 


Vq W — WqV 

kA 23 


2kA 23 wo 
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Figure 3. Minimal Squares 


A 23 = wv' — vw' = 


(wqv — vqw ) 2 + k 2 w 2 A 2 


23 


kA 23 w 0 

This means that (11311 is automatically satisfied with this choice of (v',w'). In fact, we have shown 
that, for every integer values of (v,w), such that 


(14) (w 0 v - v 0 w) 2 + k 2 w 2 A 23 = (kw '(/)A | 3 

taking (v', w') as above, we automatically get a square determined by vectors OP(v, w), P(v',w') € 
Q 4 , which may not have integer coordinates. In order for (|14l) to have rational solutions, in v and 
w, we need to have kw q£' a positive integer which is a product of primes of the form 4s + 1, and all 
the other primes, in its prime factorization, have even exponents. Clearly, there exists a smallest 
positive integer £' with this property. ■ 

As we have seen in Proposition 2.2 am we can similarly show the existence in every plane 
of equations (HOD , of a square that is minimal , in the sense that the side-lengths are the smallest 
possible. Then, a similar result to Proposition 2.2 m) takes place: every square in the same 
plane can be written in terms of a square that is minimal (same formulae as in 2 dimensions). For 
a proof we just invite the reader to study Figure [3] (a proof “without words” exercise). 

As in dimension 3, we also have parametric formulae for squares, which are more general than 
the ones at the beginning of the subsection. The two vectors are 
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u = ±(—ta — zb + yc — xd , —tb + za — yd — xc , —tc — zd — ya + xb, —td + zc + yb + xa ) 

(15) 

x = ±(ax — by — cz — dt , ay + bx + ct — dz, az — bt + cx + dy, at + bz — cy + dx) 
and one can check that u • v = 0 and |u | 2 = |x | 2 = (x 2 + y 2 + z 2 + f 2 )(a 2 + b 2 + c 2 + d 2 ). The two 
relations as in ([5]) can be simplified by (x 2 + y 2 + z 2 +f 2 ) 2 and (a 2 + b 2 + c 2 + d 2 ) 2 respectively and 
reduced to 

(x 2 + y 2 + z 2 + t 2 ) 2 = (x 2 + y 2 - z 2 - t 2 ) 2 + ( 2 xf + 2 yz ) 2 + (2 ty — 2 xz) 2 , 

(16) 

(a 2 + 6 2 + c 2 + d 2 ) 2 = (a 2 + b 2 - c 2 - d 2 ) 2 + ( 2 ad - 26c) 2 + ( 2 ac + 26d) 2 . 

These relations can be reduced even farther depending upon what q± and are. A simple 
observation here is that the simplified relations (usd depend only on the plane containing the 
square and every other square contained in the same plane has the same simplified relations. This 
implies that the minimal square contained in the plane has the sides at most y/k\k 2 , where the 
k 2 = A 2 + B 2 + Cf are the primitive versions of (fISl) . 

Example 3. We start with two equalities as in (fT6l) : 9 2 = 7 2 + 4 2 + 4 2 and 15 2 = ll 2 + 10 2 + 2 2 . 
This implies that 45 2 = 33 2 + 30 2 + 6 2 = 35 2 + 20 2 + 20 2 . As in Theorem 12.81 we can take a = 33, 
a 1 = 35, 6 = 30, b' = 20, c = 6 and d = 20. Then A12 = 1, A34 = 34, A13 = 5, A24 = 25, A14 = 13, 
and A23 = 7, and the equations of the plane are 34x + 25rc + 7t = 0 and 

7u + 5v + w = 0. The system can be solved easily in terms of u and v: w = —7 u — 5v and 
t = 25 u + 13v. Then two vectors that generate the lattice in this plane are (1,0,—7, 25) and 
(0,1, —5,13). Then the “volume” of the fundamental domain is 45. The quadratic form is Q(u, v ) = 
675u 2 + 720 uv + 195x 2 or Q(u,v) = 15(45w 2 + 48m; + 13x 2 ). So, the equation Q(u,v) = 45£ is 
equivalent to (13x + 24 u) 2 + 9u 2 = 3(13)A Then it is clear that v has to be a multiple of 3 and 
if we take t = 3, v = 3 and u = —2 we obtain a solution. This gives u = (—2, 3, —1, —11), and 
v = 6 and u = —3 gives v = (—3, 6 ,—9,3) = 3(—1,2,—3,1) (which is indeed a square). This 
means that the Ehrhart polynomial of this square is E\j(t) = 3 1 2 + 4t + 1 = (t + l)(3t + 1), t € N. 
Taking the parameterizations as in Theorem 12.41 for the two equalities, 9 2 = 7 2 + 4 2 + 4 2 and 
15 2 = ll 2 + 10 2 + 2 2 , we observe that this square is covered by the parametrization (11511 by taking 
a = x = z = l,b = d = t = 2, c = 0, and y = 3. It is surprising that such a square with a big size 
for its side lengths has no lattice points in the interior. 

In order to understand better what is happening we will reformulate everything in terms of 
quaternions. 

We remind the reader that the Hamilton quaternion algebra over the real numbers, denoted by 
1HI(M), is the associative unitary algebra given by the requirements: 

(I) IHI(M) is the free K-module over the symbols i , j, and k, with 1 the multiplicative unit; 
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(II) i 2 = j 2 = k 2 = —1, ij = —ji = k, jk = —kj = i and ki = — ik = j. 

If q = x + yi + zj + tk € H(R) the conjugate of q is q = x — yi — zj — tk and the norm of q is 

N(q) = x 2 + y 2 + z 2 + t 2 . Some standard notation is then naturally appearing: Re(q) = x and 

Im(q) = yi + zj + tk. 

By H(Z) we denote the subset of quaternions whose components are all integers. We imbed Z 4 
into H(Z) in the natural way: (x, y,z,t) x + yi + zj + tk. Also, we will think of R 3 imbedded in 
H(R) in a, more or less, natural way ( y,z,t ) yi + zj + tk\ in other words, M 3 is the hyperplane 
Re(q ) = 0. 

It is known that this norm is multiplicative, i.e. N(qiq 2 ) = N(qi)N(q 2 ), and 7 jTq 2 = 02 Ql- 
For the important results about the arithmetic of H(Z), we recommend the reader the recent 
treatment in [4]. Using the same terminology as in [4], a quaternion q is called odd, if N(q) is an 
odd number. 

The parametrization in Theorem 12.41 is basically equivalent to 

(17) n\i + rt 2 .) + n^k = q(e)q, where q = x + yi + zj + tk, and e € {i,j, k}. 

In case e = k, n\ = 2 (xz + yt ), n 2 = 2(zt — xy), and = x 2 — y 2 — z 2 + t 2 . 

The parametrization (1151) is equivalent to 


(18) u = gi£ 2(/2 and v = qie^, where q\ = x + yi + zj + tk, q 2 = a + bi + cj + dk, 

with 62 and 63 fixed in {i,j, k}. In what follows we are going to take 62 = j and 63 = k. This square 
is in a plane defined by equations similar to (1101) obtained from the two Pythagorean quadruples 
defined by q\ and qo as in (1171) using e = ei such that {ei, 62 , 63 } = {i,j, k}. 

Observation 1 : If we substitute q = q{a + j3k) into (fT71) for e = k, we see that because ( a+/3k)k = 
k(a + /3k), we have 

nii + n 2 j + n 3 k = q{k)(a + /3k)(a - /3k)q = ( a 2 + (3 2 )q(k)q = ( a 2 + f3 2 ){n{i + n^j + n^k). 

This implies that rii = ni(a 2 + f3 2 ), i = 1,2,3, which means that the primitive solutions of the 
equation Q have to arise from quaternions q which have no right factors of the form a + j3k. We 
obtain the same conclusion if q = q{\ + i) or q = q{l + j), because (1 + i)k = — j(l + i) and 
(1 + j)k = i( 1 + j). We have the following converse of this observation. 

Proposition 2.9. If in G 3 ’ we have gcd{ni,n 2 ,n^) = n> 1, i.e. if n\ = 2 (xz + yt), n 2 = 
2 (zt—xy) andnj, = x 2 — y 2 — z 2 +t 2 , are divisible byn, then q = x+yi+zj+tk = (. x'-\-y'i+z'j+t!k)rj 
where rj € {1 + i, 1 + j, a + k/3} for some integers a and f3. 
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Proof. Let us pick a prime p which divides n. Since 

n\ + n 2 = 4 [(xz + yt) 2 + (zt - xy ) 2 ] = A(x 2 + t 2 )(y 2 + z 2 ) 

is divisible by n 2 it is divisible by p. First, we assume that p > 2. Then p divides either A := x 2 + t 2 
or B := y 2 + z 2 . Since = A — B and p divides 77 . 3 , we must have, in fact, A and B both divisible 
by p. If p = 2, then A and B are either both even or both odd. If A and B are both even, then 
we have the same conclusion as in the case p > 2. In the case A and B are both odd, we have, lets 
say x = 2x', y = 2y', z = 2z' + 1, t = 2t! + 1. This implies that 

q = x+yi+zj+tk = 2 x'+2y'i+2z' j+2t'k+j+k = {x'+y'i+z' j+t'k)(\ — i)(l+i)+k(l+i) = q(l+i) 

which proves our claim. Similar conclusion can be drawn if x and z, t and y, or t and z are even. 

If p is a prime of the form 4 k + 3, then automatically x = px ', z = pz', t = pt !, and y = py', 
which shows that q = x + yi + zj + tk = (x' + y'i + z'j + t'k)p and the conclusion of our proposition 
follows, with a = p and (3 = 0. 

Finally, if p is a prime of the form Ak + 1, then p = a 2 + /3 2 . Because A = (x + tk)(x — tk) and 
y = a+/3k is a Gaussian prime integer, it divides x+tk or x — tk. Without loss of generality we may 
assume that y divides x + tk. If this is not the case we continue with y. Next we observe that since 
p divides xz + yt and zt — xy , then p divides (x + tk)(y + zk ) = (xy — tz ) + (ty + xz)k. This implies 
that if y does not divides x + tz, then rj must divide y + zk. Since y divides B = (y + zk)(y — zk), 
then y divides either y + zk or y — zk. Let us assume first that y divides y — zk. Then we can write 

q = x + yi + zj + tk = x + tk + i(y — zk) = (x + t'k)y + i(y + zk)y = qy. 

If y divides y + zk, then rj divides y — zk. If rj divides z + tz then we proceed as above, with a 
small change: 

q = x + yi + zj + tk = x + tk + i(y — zk) = (x' + t'k)y + i(y' + z'k)y = qy. 

If rj does not divides x + tz, then we have shown above that in this case, rj must divide y + zk, 
which is the same thing as y dividing y — zk. This puts us in the same position as above. ■ 

Observation 2: If we substitute q\ = qi(a + /3i) into (fTKj) . we see that because (a + /3i)k = ak — f3j 
and (a + f3i)j = aj + (3k, we have u = au — f3v and v = f3u + av. This is saying that the twin 
pair u and v is in the same plane as u and v. A similar statement can be obtained if we substitute 
<?2 = < 12(01 + (3i) into CED- 

Let us recall, for the convenience of the reader, the statement of Lemma 2 . 6.5 in [ 4 ], with the 
only difference that the important factors appear on the right (so, this new statement follows by 
conjugation and a symmetry type transformation i —>• —i, j —>• — j, and k —>• —k). 
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Lemma 2 . 10 . Every integer quaternion q £ H(Z) has unique factorization 
(19) q = 2 % V , 

where q € H(Z) is odd, rj q € { 1,1 + i, 1 + j, 1 + k, (1 + j)( 1 + i ), (1 — fc)(l + i)} and for some 
non-negative l q € Z. 

This new statement follows by conjugation and a symmetry type transformation i —>• —i, j —j , 

and k —>• —A; from the original statement in [4|. 

Theorem 2 . 11 . (i) We assume that q\ and q 2 m f/ie parametrization U8\) represented as in 
m , are not right-divisible by quaternions of the form it = a + ( 3 i, | 7 r| > 1 , then the square in the 
parametrization m is minimal. 

(ii) The parametrization m represents all the integer squares. 

Proof, (i) We observe first that under our hypothesis, i qi = i q2 = 0 and rj qi , r) q2 £ {1+j, 1+k} 
in (1191) . By way of contradiction, if the construction of the square S = {u, u} in (1181) is not minimal, 
there exists an integer square So = {no, no} in the same plane in such a way S can be obtain from 
Sq by a simple transformation which involves two integer parameters a and /3: u = auo + /3v 0 
and v = —/3uo + avo- We may assume that p = a 2 + f3 2 is a prime, otherwise we redefine the 
twin pair Sq in such a way this condition is satisfied. Let us first assume that p > 2 . Since q\ 
is not right-divisible by rj = a + (3i (which is a prime in H(Z), i] odd quaternion), we conclude 
that gcd(< 7 i, rj) r = 1 (the right greatest-common divisor as defined in [3]). By Theorem 2.6.6 in [3j, 
which is the equivalent of Bezout’s relation in H(Z), we can find 7 and 6 in H(Z[^]) such that 

7 <Zi + Sfj = 1. 

If p = 2, then \u\ 2 = \v\ 2 = N(qi)N(q2) = 2|«o | 2 and so 2 divides N(q\) or N(q 2 ). Using 
Lemma f 2.101 we conclude that q\ is of the form q\ (1 + k) or q\ (1 + j) with q[ odd. So we can still 
obtain a Bezout’s relation as above, with q[ and rj, which can be easily transformed into one for q\ 
and rj. 

Solving for uq, the above system, in terms of u and v, we have 

u 0 = (1 /p){au - / 3 v) = (l/p)qi(a + i/ 3 )kq2. 

Multiplying on the right with (1 /p)(a + if3)kq2 the above Bezout’s relation, we obtain 

yito + dkq^ = (1 /p){a + i( 3 )kq o. 

From this relation we see that q2 = < 72(0 — i/ 5 ) for some <72 € H(Z[-(]). Hence 2 k q2 = ^ 2 ( a ~ */^) 
for some q' 2 £ 1 HI(Z). Using Lemma [2.101 again, we conclude that q2 is right divisible by a — i /3 which 
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is in contradiction with our hypothesis. Therefore, it must be true that the S = {u, v} in (1181) is 
minimal. 

(ii) Let us start with an arbitrary square, So- We then use Theorem 12.71 to find the reduced 
equations of the plane that it is contained in, as those given by m- Then we can apply the 
Theorem 12.41 twice, and finally use the two sets of parameters (8 in total) to construct a square <Si, 
with formulae (1151) . The quaternions involved have the property in part (i) (both odd), and so the 
square 5i is minimal. Then can be written in terms of 5i using two numbers a and (3. Using 
the substitutions described before the statement the theorem, we can we can write 5o as in m 
by multiplying either of the quaternions by a + / 3 i. ■ 

Theorem 2 . 12 . Assuming we have two odd integer quaternions q\ and c/2 as in (ii) in Theo- 
rem \2. lA then the fundamental domain of the integer lattice containing the square {u,v} in U8\) . 
has “volume” equal to V = lcm(N(q\)), Nfa))- The Ehrhart polynomial of the square U8\) is 

En(t ) = gcd(JV(gi), N(q 2 ))t 2 + (D 1 + D 2 )t + 1, 
where D\ = gcd(ui, u 2 , U3, U4) and D2 = gcd(ui, v 2 , V3, V4). 

PROOF. This follows from Theorem 12.71 (part (ii)), Proposition 12.91 and the observation that 
lcm(N(qi),N(q 2 )) 2 = a 2 + b 2 + c 2 = a ' 2 + b ' 2 + c 12 for some integers a, b, c, a', b', and d with 
gcd(a, 6, c, a', 6', d) = 1 . ■ 

Observation 3 : To get back in dimension 3, we identified M 3 with the subspace of quaternions q 
for which Re{q ) = 0. We have seen that only one Pythagorean quadruple equation is necessary: 
A 13 = A 12 = 0 which implies A 14 = 0. To accomplish this, we can take q\ = c /2 = q where q is 
odd. This implies that the most general square in dimension 3 is of the form Q combined with 
the usual variations from dimension 2 in the respective plane: u = au — /3v and v = f3u + av with 
{u, v} as in Q. 

There seems to be plenty of numerical evidence that the sequence of almost perfect squares in 
dimension four is the set of all non-negative integers. Let us denote the set of all such numbers by 
APS4. We include some partial result in this direction. 

Theorem 2 . 13 . (i) The set APS4 contains all odd numbers. 

(ii) The set APS4 contains all even numbers of the form p — 1 with p a prime p > 11. 

Proof, (i) Let us take an odd k, k > 1, and a primitive Pythagorean quadruple representation 
k 2 = a 2 + b 2 + c 2 . Without loss of generality we may assume that a is odd and both b and c are 
even. If we multiply this equality by 2, we obtain 2 k 2 = 2a 2 + 2 b 2 + 2c 2 = 2c 2 + (a + b ) 2 + (a — b) 2 . 
Hence we can build the square u = ( k , k, 0,0) and v = (c, —c, a + b, a — b). It is clear that D\ = k 




EHRHART POLYNOMIAL FOR LATTICE SQUARES, CUBES AND HYPERCUBES 


19 


and D 2 = 1 (a — b and a + b are odd) as defined in Theorem 12. 121 Since the sides have length kV2 
and the two equation defining the plane containing {it, v } are k 2 = (—c ) 2 + b 2 + a 2 = (—c ) 2 + a 2 + b 2 , 
we conclude by Theorem 12. 121 that the Ehrhart polynomial is E UjV (t) = 2 kt 2 + (k + 1 )t + 1. Hence 
E u ,v{— 1) = k. This proves the first statement. 

(ii) If p is a prime p > 11, then we can find two distinct primitive Pythagorean quadruple 
representations p 2 = a 2 + b 2 + c 2 = a ' 2 + b ' 2 + d 2 and construct a square as in (fT 8 l) , whose Ehrhart 
polynomial is Eft) = pt 2 + (D 1 + D 2 )t + 1, where Di divide p. It is not possible to have = p 
since this puts us in the position where one of the vectors defining the square, say it, is of the 
form u = ( p , 0, 0,0). This leads essentially to only one equation which defines the plane containing 
{it, v}, but we have two different equations that define the plane. Hence, since E{— 1) = p — 1 the 
claim in (ii) follows. ■ 

Numerical evidence shows that given an odd k > 9, one can find a square whose Ehrhart 
polynomial is Eft) = kt 2 + 2 t+l. For instance, for k = 2015 we found it = (—836,584, —1592, —697) 
and v = (—506,1414, 203,1328) that does the job. 

2.4. Cubes. If a general lattice cube in K 3 is given by the orthogonal matrix 

ai bi c\ 

b 2 c 2 , 

03 b% C 3 

with ai, bi and c,; integers satisfying OjOj + bibj + CiCj = 5ij £ 2 for all i, j in {1,2,3}), we define 
di := gcd(aj, q). It is clear that the di are divisors of £. Assuming that the matrix in (1201) 
is irreducible (no smaller £ can be used, which implies l odd), then we have shown in |18j the 
following. 

Theorem 2 . 14 . Given a cube Cg constructed from a matrix as in M), its Ehrhart polynomial 
is given by 

£ 3 t 3 + £(d\ + d 2 + dz)t 2 + ( d\ + d 2 + d^)t + 1 or 
( 21 ) L(Cg,t) = ,t€ N. 

[it + 1 )[£ 2 t~ + (di + (^2 + ^3 — £)t + 1] 

Similar formula takes place for the cubes in dimension four. In general a square can be completed 

in various ways to a cube and even to a hypercube. Let us assume that the cube is given by the 

first three vectors in the rows of the following matrix 

<21 61 Cl d\ 

a 2 b 2 c 2 d 2 
03 b 3 c 3 d 3 ’ 

O 4 64 C 4 C ?4 


( 20 ) 


C ‘ = l 


(22) 










20 


EUGEN J. IONASCU 


with Oj, bi, Ci and di integers aiaj + bibj+CiCj + didj = S t j£ for all i, j in {1,2, 3,4}). In this case, we 
define Di := gcd(aj, 6 *, a, di) and Qj be the greatest common divisor of all the 2 -by -2 determinants 
of the matrix 

CLi bi Ci di 
a j bj Cj dj 

for every i and j with 1 < i < j < 3. 

Theorem 2.15. Given a cube C £ constructed from a matrix as in if 221 ) . its Ehrhart polynomial 
is given by 

(23) L(C e , t) = IDyt 3 + At 2 + A't + 1 , t € N, 

where A := Ci2 + Ci3 + C23 and A' := Di + D 2 + D 3 . 

Proof. For the first coefficient the volume of the cube is l\ft and the volume of the fundamen¬ 
tal domain of the lattice containing it is equal to-^/a 2 + b\ + c\ + d^/D^ = jy^. For the second 
coefficient, the area of each of the six faces is t. Using the Theorem 12.111 to find the “volume” of 
the fundamental domain of the lattice containing the faces determined by rows i and j, we obtain 
£ / S t j , which implies the stated value of A. Finally, for the coefficient A', we proceed like in the 
proof of Theorem 12.141 ([18]). The main idea is essentially based on the fact that the topological 
equivalent in our cube to [ 0 , l ) 3 (containing k lattice points) has the fundamental domain property: 
the dilation with a factor t contains exactly t 3 n lattice points. This implies that k = ID 4 and so 


£D± — 1 + (ID^ — A + A / — 1) + ^ {Cij ~ (Di + Dj) + 1) + ^ ^(Dj ~~ !)■ 




Solving for A' we obtain the stated value in the theorem. 


2.5. Hypercubes. Hypercubes in dimension four can be constructed using quaternion techniques 
described earlier. A few orthogonal matrices in dimension 4, together with the Ehrhart polynomial 
associated to the corresponding hypercube, are included next: 


■ 1 

1 

0 

0 


1 

1 

1 

0 


1 

1 

1 

-1 ~ 

1 

-1 

0 

0 

1 

-1 

1 

0 

1 

1 

-1 

1 

1 

1 

0 

0 

1 

1 

’ 71 

0 

-1 

1 

1 

’ 2 

1 

-1 

1 

1 

0 

0 

0 

-1 

-1 

0 

1 

-1 


1 

1 

-1 

1 


1 

72 


4 1 4 + 8t 3 + 8t 2 + At + 1 
= (2f 2 + 2t + l) 2 


9 f 4 + 12t 3 + 6t 2 + 4t + 1 
= (t + l)(3t + l)(3t 2 + 1) 


16f 4 + 16t 3 + 12t 2 + At + 1 
= (1 + 2 1 + 4 1 2 ) 2 
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' 2 

1 

0 

0 


' 2 

1 

1 

0 


'21 1 1 

1 

-2 

0 

0 

1 

1 

-2 

0 

1 

1 

1 - 2-1 1 

0 

0 

2 

1 

’7! 

1 

0 

-2 

-1 

’7f 

11 - 2-1 

0 

0 

1 

-2 


0 

1 

-1 

2 


1 - 11-2 


25i 4 + 20 1 3 + 14 1 2 +41+1 36t 4 +24i 3 +8t 2 +4i+l 49++28t 3 +6t 2 +4t+l 

= (1 + 21 + 51 2 ) 2 



3 

0 

0 

0 


2 

2 

1 

0 


2 

21 1 

1 

0 

2 

2 

1 

1 

2 

-2 

0 

1 

1 

2 

-2 -1 1 

3 

0 

2 -1 

-2 

’ 3 

1 

0 

-2 

-2 

’ \/Ib 

1 

1 -2 -2 


0 

1 -2 

2 


0 

1 

-2 

2 

1 

-1 2 -2 


811 4 + 541 3 + 181 2 + 61 + 1 81t 4 +36t 3 +6t 2 +4t+l 100t 4 +40i 3 +16t 2 +4t+l 

= (31+1) 2 (91 2 + 1) 



3 

-1 

0 

0 


3 

1 

1 

0 


2 

2 2 1 

1 

1 

3 

0 

0 

1 

1 

-3 

0 

1 

, and ,_ 

\/l3 

2 

-2 1 -2 

vTo 

0 

0 

3 

1 

’Vn 

1 

0 

-3 

-1 

2 

-1 -2 2 


0 

0 

1 

-3 


0 

1 

-1 

3 


1 

2 -2 -2 


1001 4 + 401 3 + 241 2 + 61 + 1 121t 4 +44i 3 +6t 2 +4t+l 169i 4 +53t 3 +6t 2 +4i+l 

= (101 2 + 21 + l) 2 


We see that several of these examples are cross polytopes (squares in two dimensions) and indeed 
the resulting polynomial is the product of the smaller degree polynomials involved. Let us assume 
that in general we have 


(24) 



a\ b\ ci d\ 

GQ &2 C-2 d2 

«3 ^3 C3 G?3 

04 64 C4 C?4 


where <+aj + bibj + CiCj + didj = 8^1 for i. j € {1,2, 3,4}. We define Di = gcd(aj, bi, Ci, di) for 
i € {1,2,3,4} and assume that H£ is irreducible, i.e., gcd(Z?i, D 2 , D 3 , D 4 ) = 1. The Ehrhart 
polynomial associated to the hypercube constructed in a natural way using the orthogonal matrix 


Hi is then 


(25) Eh (£•)(£) = + ait 3 + Q! 2 1 2 + a$t + 1. 

Proposition 2.16. With the notation introduced earlier a \ = i{D 1 + D 2 + H3 + D 4 ). 

Proof. Let us look at a 3-dimensional face, say generated by first three rows of (1221) . Its volume 
is iVi and the volume of the lattice containing it is equal to\/ a\ + 64 + c\ + Since 

there are eight such faces the general theory gives the value claimed. ■ 
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Theorem 2.17. We have 


(26) «2 = as + 5 — A, 

where 6 = Yj\<i<j <4 Cij an d A = D\ + D 2 + D 3 + Z) 4 . 

Proof. We are going to use the same property of the part of hypercube that is topologically 
equivalent to [0, l) 4 , as in the proof of Theorem 12.151 The balancing of points in each of the 
corresponding sets of the partition 

[o, i ) 4 = (o, o, o, o) u (o, i ) 4 J (O ’ 1 ) 3 U (O ’ 1 ) 2 U (O’ 1 ) 

all four 3 D cubes all six 2D faces all four ID sides 

gives (using Theorem 12.151 ) 

= 1 + ( J? 2 — oil -\- 0.2 — oi 3 + 1) + {IA — 25 + 4A — 4) + ^ ^ (Cij ~ {Di + Dj ) + 1) + A — 4. 

From this, one can easily derive (1261) . ■ 

Looking at the examples we have so far, we notice that a 3 = A, and so the Ehrhart polynomial 
of an hypercube would result into the follows simple form 


(27) E H{e) (t) = ft 4 + £A t 3 + 5t 2 + At + 1. 

This is indeed the case since it follows from Theorem 9.9 in [2]. Of course, one may want to 
generalize all these concepts to dimensions bigger than four. 
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